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000000 TITTTTITTT SSSSSSSS PPPPPPPP 000000 wat Ww HH HH HH HH 
000000 TITTTITITT SSSSSSSS PPPPPPPP 000000 wi WwW HH HH HH HH 

00 00 TT $$ PP PP 00 00 WW WW HH HH HH HH 

00 00 TT SS PP PP 00 00 WW Ww HH HH HH HH 

00 00 TT SS PP PP 00 00 WW WW HH HH HH HH 

00 00 TT SS PP PP 00 HH HH HH HH 

00 00 TT SSSSSS$ PPPPPPPP 00 00 WW WW HHHHHHHHHH HHHHHHHHHH 
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00 00 TT SS PP 00 OO WW WW WwW HH HH HH HH 

00 00 TT SS PP 00 00 Ww Ww WwW HH HH HH HH 

00 00 TT SS PP 00 00 WWWW WwWWwW HH HH HH HH cece 

00 00 TT SS PP 00 00 WWWW WwWWwW HH HH HH HH cece 
000000 TT SSSSSSSS PP 000000 Wt Ww HH HH HH HH cece 
000000 TT SSSSSSSS PP 000000 wid Ww HH HH HH HH eevee 

LL III! SSSSSSSS 

LL III] SSSSSSSS 

LL I] S$ 

iL I] SS 

LL I] SS 

LL I] SS 

LL I] SSSSSS 

LL I] SSSSSS 

LL Il $$ 

LL I] SS 

LL I] SS 

LL I] SS 

LLLLELLLLL III] SSSSSSSS 

LLLLLLLLLL HII! SSSSSSSS 


oO 
Oo 
= 
=< 
<= 
< 


} 


abl 


TSSPOWHH 
e of contents 
(2) 4 
(¢) 3 
(4) 2 


C 
= REAL®16 ** REAL*16 power routine 


ORY ; Detailed current edit history 
ARATION 
POWHH_RS = H floating ** H_floating 


16-SEP-1984 02:00:37 VAX/VMS Macro v04-00 


Page 


0 


| D 2 | 
TS$POWHH ~ REAL*16 #* REAL#16 power routin 16=SEP=1984 02:00:37 yAX/VMS v04-00 Page 1 
806 PP apes nSEP=198¢ 14:58:21 EMTMRTL.secdorsPoumN.maR:1 “29° cf) 
ne 
IT: 


¢ 
SRC 
; 9 SSP + ae - pee ** REAL*16 power dong * 


TS 
2-006/ File: OTSPOWHH.MAR E JCW2006 


oooo 


OOoOoooooooooooooooo°oo 


CGOOCOCCOCOSCSCOSOSSSCOCCOSOSCOCOSCOSOOOOOOOOCOOOOOOOOOO 


Dahebhh bh RRRRRRRRRRRARSSAALESALERSASASEL SESSA EAR R ARES ER SESE RRR REESE EEE SEES 


COPYRIGHT (c) 1978, 1980, 1982, 1984 

DIGITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 

THIS SOFTWARE Is FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
N ACCORDANCE WITH THE TERMS SUCH LICENSE AND WITH THE 


T OT cE 0 
THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE A iy ee IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
ooePokat ion. NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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: FACILITY: Compiled code support Library 
ABSTRACT: 
H_floating base to H_floating power 
VERSION: 2 


AUTHOR: 
Bob Hanek, 9-Mar-83; Version 2 
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TSSPOWHH - REAL*16 ** REAL*16 power routine 16-SEP-1984 02:00:37 \VAX/VMS Macro v04-00 P 
O86 wIstoRy ; Detailed current edit history $78 eats ri 94:98:34 MTHRTL.SRCJOTSPOWHH.MAR; 1 — 


ang . =fABLE ey to avoid Linker problems with .ADDRESS Tor data. 
-FRay> : 
$-003 Change rene ining reference of INDEX(Rx) to INDEXCRxJ. LEB 29-May-1983 
-004 Add in # signs to_avoid Linker errors regarding non-relocatable 
references. LEB 30-May-198 
$-00? Correct use of quadword index from INDEX. SBL 31-May-1983 
-006 Corrected a bug involving a SYS_F_FLTOVF_F error dur "9 3 
MULH2 R4, t2(SP). Code was added fo see 7f a MTH overflow message 


44 ~SBTTL HISTORY ; Detailed current edit history 
; Edit history for OTSSPOWHH 
: 3-001 Implemented new algorithm. RNH 18-Mar-83 
; 2-002 Chan ; Taatetaes Of A1_TABLE(Rx) and A2_TABLE(Rx) to A1_TABLECRx] 
38 
000 
0000 
0000 
4464 or a zero should be returned. JCW 19-Jan-1984 
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* REAL*16 power routine 16-SEP-1984 00:37 VAX/VMS Macro Vv04-00 P 
§ . greee 188s 7:5R531 EAtAKTe. eaesorseoung.mar:1 2% (3 
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) 


; INCLUDE FILES: 


: EXTERNAL SYMBOLS: 


| 
| 
DSA GBL 
eEXTRN MTHSSSIGNAL ; Math error routine 
~-EXTRN MTHSK_FLOOVEMAT Floating point overflow code 
Tih. MTHSK_FLOUNDMAT Floating point underflow code 


RN MTHSK_UNDEXP Undefined exponentiation code 


: MACROS: 
: SSFDEF :; Define stack frame symbols 
: EQUATED SYMBOLS: 
: base = 3; offset to base (by value) 
exp = 20 3; offset to exponent (by value) 
ACMASK = “M<R4, R5, Rb, R7, RB, ROD 
; register saving mask 
t2 = 4 ; Offsets from stack pointer | 
t4 = 20 ; to temporary storage locations. 
t6 = 36 : Each location is 16 bytes 


OOOO 00000009 090909 0 O09 NI SINISE NINO A OAOAOOAO 
FW 9 OONA UE WN 9 OONAUE WN OC OONOAUE UT 
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TSSPOWHH REAL*16 ** REAL*16 power routine 1 P=1984 7 VAX/VMS Rs cro V04-00 Page 4 | 
Br 806 DECLARATIONS : rae: Pa} 984 tf 98: 3 MTHRTL.S REI TSPOWHH.MAR; 1 7 (4) | 
: 
6 3) ; PSECT DECLARATIONS: 
000000 39 -PSECT _OTSSCODE PJ.C,SHR,QUAD,EXE,NOWRT 
88 199 3 program section for OTS$ code | 
444 199 ; CONSTANTS: 
00 105; 
on8 104 
00 193 ; 
0000 1 $ 3; The INDEX table gives the offset | B gyecyerds? to the a pce riate 
8898 18 ; entries in A1_TABLE and A2 “Tein E: Entry 1 of the table | 
108 ; is a special @ncoding that” laeaded to access the octawords immediately 
0000 109 ; BEFORE the Al and A2 tables. 
$33 e 
06 04 04 04 02 02 FF 00 $008 M6 INDEX: .BYTE “X00, “XFF, “x02, xOee “x04, “X04, “X04, *X06 | 
OA OA OA 08 08 08 06 06 0008 = 11 -BYTE *X06, “X06, “X08, “x08, “X08, “*XOA, “XOA, “*XOA 
10 OE OE 13 of Of 0¢ OC 0010 114 -BYTE ‘*XOC, “XOC, “xOC, nt “KOE, “XOE, “XOE, “x10 
14 14 14 12 12 12 10 10 0018 115 BYTE “x10, *X10, *X12, *X12, *X12, *X14, *K14, °X14 
18 18 18 18 16 16 16 14 0020 116 -BYTE “X14, “X16, *X16, “X16, “X18, *X18, “X18, *x18 
1€ 1€ 1C 1C 1€ 1A 1A 1A 0028 117 «BYTE “*X1A, “X1A, “X1A, *X1C, “K1C, *X1C, “X1C, “XIE 
22 20 20 20 20 a5 45 45 0030 118 BYTE atd5- aXt5. *K1E, “X20, *X20, *X20, “X20, *Xx22 
26 24 26 26 26 0038 119 -BYTE “*X22, “X22, *X22, *KX24, “X24, *KX24, “X24, *X26 
28 28 28 28 28 26 26 26 0040 120 -BYTE *X26, “X26, “X26, *X28, “X28, *X28, “X28, *X28 
$5 C 2C 2C 2A 2A 2A 2A 0048 121 BYTE “X2A, “K2A, “K2A, “KX2A, “X2C, “X2C, *KX2C, *X2C 
0 0 0 5 5 3 f f 0050 1 ¢ -BYTE | es “x - “x 5. “x 5. “x # = 0, “x 0, “Xx 0 
34 34 32 0058 12 «BYTE 2350, “R50, “Se. *XSe, “Sa, *NSe. *XS4, "854 
36 36 36 36 36 34 34 34 0060 124 -BYTE “X34, “X34, “X34, “X36, “X36, *X36, *X36, *X36 
3A 3A 3 38 38 38 3 8 0068 125 -BYTE “*X38, “X38, “X38, “X38, “X38, “X3A, “X3A, “X3A 
3c 3c 3c 3C 3C 3A 3A 3A 0070 126 -BYTE SRSA, “ASA, “ASA, ASE, MRSC. TUS, MAS. TES 
40 40 40 3€ 3€ 3E 3E SE it \st BYTE “X3E, “X3E, “X3E, “X3E, “X3E, “X40, “x40, *x40 
0080 \¢2 ALIGN QUAD 
0080 8130 
0080 131; 
0080 13 3; For k = 0, 32, the k-th entry of A1_TABLE is value of 2*(k/32) 
0080 133; rounded to 1 its" fraction bits and the k-th ehtry of A2_ TALE 18, the value 
0080 134 ; of 2(k/32) = Al_TABLE(k) rounded to 113 bits. For k = Al ABLE gives 
0080 135 : the value 2*(1/84) rounded to 113 bits and A? _TABLE give ‘gat t 760) 
Boe : § ; A1_TABLE(-1) rounded to 113 bits 
080 138° 
TA2ZACAGF F7CAOEE6 7806A3E7 02C94001 34 172 i ns “X7A2ACAGF F 7CAOEE67806A3E 702094001 
0900000 0 04001 0090 141 ~ .OCTA *x000000000000000 000 001 
GANBYIB 85 avetey ae 544 B3570081 039 sie BAO 166 -OCTA a8 ‘® sare 9 seas E3837 Ost 19998¢001 
8B92F 62 25631 444 0B0 14 -OCTA “*X4 F0B554001 
8DAED9 OOF OE ee 2 be 01 0CO 144 -OCTA *X3 BDAED9BRr {0 27620988 561i4 04001 
Mh ratat o EE 9 84001 0D0 145 -OCTA ae aS fF cea o6p ; €7172B4001 
oF OR <a?" SBBOA 80 ah e? 33 om 1 BED 196 -OCTA anerarete 5880A780 38 1610484 4 
5CB6B1C1 1FA Z3AbE? 28 eee 01 F 14 -OCTA “*X5CB UCI IFADS gc pare 874001 
rarret yt sb 12C 443434 1 010 168 -OCTA “*X4AA4F5A25D1 He 44473: 94001 
2686465 D5 1B if A ty 1 11 14 -OCTA “nis oe? 3 AS HH 5 A 9 F4001 
4 et 455 AA9 A4001 190 139 -OCTA “*X4 F4 gh 9¢737 A4001 
077435E13 412 e338 $38 ett DEASOO1 0130 151 -OCTA *XD7745E1341 Et 426 4 DEA4001 
O1AOO9DF 36F4D031 18928606 44£04001 0140 152 -OCTA “*XO1AOO9DF 36F 40031 33 860644E04001 


16° 


2 


power routine 


= REAL*16 ** REAL*16 


DECLARATIONS 


TSSPOWHH 


2-606 


SEp=1986 i:B:21 ERIMRTL.SRCSOTSPOUMN.mAR:1 "28° 


eR NU / Oou Owuus Ou Ou OVu wu We VOWWEuenunuuuu Or 
SOOOOOOOCOCOOGOOOOOOOOOOSoSSoO @® ©0000 60.60 60 00 00 00 00 60 60 G0 60 50 00 00 00 00.00 G0 00 00 G0 00000060 60000000 
SOooOooooooooooocoooeooo uu Output uw_uuwuuwuwuuw a Oe Oe te A Oe A 
> > > > > > > > >>> >> Oe) O17 OP) OO P90 9 OP) 1) OD PEE OPEN 
QMO. OF. -DuUMOOSROUK Oro @D OWEMON se OF DOM“ MODMNO LOOT OOM NO OOu 
Eh A pee te etter epee tae S64 wm OwWO LD OMNM & FOP PMID & WOU. OOM UNOS Ou 
Om ave aVMOMNUw- @OMiu NO oo SLUM TA OM OF COP OPN~T ND OPI Wu Lu COCD00 
TUNOOP-F. DODOO Mv VAawuod SO o~unwwh we Ow OMe & oR We OU UE NOOMU a O 
MOTI wus VO MN DusurDoo a OF ENN Ee OW WI MOO Ow OMVLOACOe @amow 
moa MAVLMUNDAO MM OM LOTR WT OTNOODVE FE TOOuUNMMWwWOunm 

OMNADOuUlO TF DOOOuUOOQamMano < Se ene tee wm DOWOOWN £oO 
OP WUANTOU—- Mu OCODMOOru0O O&O SOMAMOAULMOM Sees ENnw et LUO LOOw TT MNo tw 
CUO VP. DW OWN ETO ODOMATO WwW Oor-OC wh uh OP & UME =P UMINAAA ODOM 
Ke FTNOu —WOMOUMWUVOTDODOMO 2a SUMO QwWOOOww CVOMA0CwWOO +s fC MOWUMOD 
MTMNUWWNOVU LUUWVaUNMO aOowTo ws OWAWe OM Ou. OO << CMNM QMO BP MOMOu tu 
SR SS BuOOrwo Pr OOaTth FT OOO UMW FVOOU Mu & — OOAWVOC OOu 
TOVUDVOV aK. TOaVveavoeaysTo a OD OA0— OM Tu OWDOKP— OOM - ODQAWNVO Ou OKO 
Qau.@0-OoOrh. ff OTVOwaewmaver.o @® Ost LOT VO LOO OO 2O Ou ST OM POO Tn 
—AMNDO oO AdMNu — WN VUMODOO ww SUM A09000 COW st — WO Wie LVOWOO OOF WOW 
LK ATF OUMOVUOTOVAOVORWAOCOO O&O SQV NET uO Om C4 OO LOWOOO we VO OWWw 
ZV OOuvaiNm@Ooad-AamtVaCDoO —-— OTNWVDDNOu ODM AMF VOR. VLOWOUMOOusrr 2 
MMO fu COON YOTTOOOM=aIM~ TO Oo OS nee eee nw Oe eae O-wow 
OOM a OF DV TF TMODOVOLCR DO CO WOUmMorWwOoudew dt oodusrver <hrMoa our VW 
MOO @rmuf-OVUCFR.OMaEODut OO — JOSr.©OSMOC Ou <OrOwovworc aw oo — 
TO WOO VUMVT EPPING TFTOOMDO WY SI FOO WIM WIT OOO O ODN. Ou OVO OD wzro<« 
VOC 004.000 M- ar. DOC OwUNo uu OMNWO OP & OOF FTUVANSO a Tooaouumoono 
WOK ea Muse Ce OP OF EEWOODOO DO KOOMF UE MOMNWOM OW O- DOVE ONVOUWOVR- VO et 
SeTNO- Mu fe Ore K VK UM O fF oa ae fe nd eee Baye A he ee 
eT TN OM OOLMOF-AODCOSO SC —COOO0OsTO-OOoorOoouUsd Voor ws VO~wrOoOo 
MOMTOO-“OOVWOUNMEDO-K- MOO << LOrK-wor-TODNDOOOCDO-— OVC Lr VOVOR—-VOOr-uO 
wet Ot etvouet coe OrOVvMooeawo ww 1'104704uwewooomunr <£OvsOu LOOM OC <OOCOMWY 
MOWOuUst OT VONONMeK SWOOO WwW "ONO & LOK TT OOM VP. OM OFM DO OMWOOW er OW 
>< DC DC OK OK OK OK OK OK OK OE OE OE OE DE COC DE KKK OKC OK OE OK OK DK OK DOK DK DK DK DK DK DK DK PE EE DK DK OK OK OK OE EK OE OE OE OE OO OK 
eeeeeee eee eee eee ee eee < pi SPE PSEC SSESEASEHEEESAAAEIETEESEES 

z 

“qq@4adddeede eee eee ee eee eet <«  wdtde de de ede ede ede dee deed eee eee eee eee eect 
el eet cael cet eel eae ee meee ee cee cee ce el ce et ee Sl el el cee 
VBVVOVOVVOVVOVOVVVVVVVOYVYVVYU Vo DEV VVVVVVV VV VV VV VOU VOU EVO VUVUUVUUYUYU 
COOOOWOWVODOOCOOOOCOCCOCOCOO AO <O0COCOCCOCOCCOCOCOCOCOoCOoCoCooOoCoOoooo°ocoeoooo 


PMTMNOPR. DOO CUM TN OR. DOOR VUM TN TUM TMAOR. DOO K— CIM TWN OR. DOO CUM TNOR. OO 
MMAMANUANOOOOOOO OOO LLL R OGD 30.60 6D 60 CODIDBAAAAAAO OA APAOOOOOOOOOO 
NNN NUNNINIONIONN 


OOOO OOCOVOOOOOCOOCOCOCOCSCOCOSOOOOSOOSSOOCOSOOCOSCOCOOCOOOOOOOSOOOOOOOOoOoOO 
WOPR.OO OV OWE OK- CUM TNOR. OO Cf I OMOOVY OWL OK UM TNOR.OG COV QWELOK MUM TNO. DOO 
el ee ee os cucucucucucucucucy ITTV IPT EAP {woth > > &> &> > > & & > > i 
OOOO OCOOOOOCOOOCOCOCOCOOCOOCOoOOoO.: SOOOCOCCOCOCOCOOC COCO OOOOOOOOOOOOOOOOO 
tek see geek sei eet eel cel cel eek onal cel el cel cee oe cee ee see el NM & OC$“OO4u 6844.44 C4 Du OVE ee We VUWWe eee ee OF 
COOOOOOCOCOCOOCSOOOCOCOCOCOSOSSS BD WOOD DOOOD GO 60 G0 60.60 cD. 
FISIIISSSSSIIIIISSIFZI OH SSMGEM GGG OM GSS GM Ge avy ww Bere ven wee ow 
ammrnn 
ANCOR. OF.  - DUMP OPP OVO LLCO - oownu. Se ee eS en Sma Roa ont aonNew ~oaouw 
eT OTFTOUK NR TNODULPPDWOTFOO MM COOwWODQMNM & TF OMPMVUDe WOW. KO OUVUMNOW Ou 
OM<aivd—- aVaEMoUMNWE- Oui“u Ho @w SOL NNN Dee BO EONS BOM Ue ee 
TNOOP-PR.DODOO EE OMOVOOWUO OC COR UMNWMNOOWU— Ou OR uu £R We OVE EMOOMUUu 
PATI; Ww VO EWN Dur LOO A OF a NnNNw eu OMe We MOD we OOr2 Oe Oe. <mow 
o mm OM LOTR WT OTNOODVEK KR TOOUMN rey) TS 
o <« SO BAN AGAE A we TF DONODOOOu <a 
~ o oO MMOur Toenrwnnw VO<oowyT-No<dtwv 
uw ow 33 OF & Ma wr. —AWPMINAAA OAD UW 
rs o a Sum OQwWwOOO wwe OCVOMOOWOO vy << DOVUMOD 
VTMMNMNOY fvuwouw s owo ww OWnOwWu ONO VUu OO << EMMM QMNOWR. oR NOMDu 
OOO du VOveomMw<eOnuDOwwo fr SOO OTR. FOOD <UMW FVOOUR Mae — OOWVOO OOu 
TOVuODVOV au TOaQVOOOaATO a O@Ma@OnK UMN OO a Oe No rt pu 
On@-Oor <4 Owvrvaw lhe ty As | [. } ard aw | < Sree aes ¢ eFroMnwes VOT TOU 
-aw BOO Au HK BAN VDoD w uM OODOO<E wor we ew Qa Nowe 
ZV-OOuvae COa0-OadtV<40 — a eRolsams rs) ewJw—rS 
SOAKNK HOOK SIMBOCOONDS B Sunoswocoeusaorwnarar SS2SRRN2S 
— . ad - “ Cs io - 
AVoSsamnerounnromzeue - SxS SRSs ee Se Sisson <<casn DMM 
won QMS te A MUA awe =) WETOO WOM We OOOO ODM Ou OYVOADYOOvO< 
SERSS= 5 Ar <P aos owen re Aneon TTVNCONGE Nee VWMAOO 
WA eM Se CUR. OF. < LWOODOO @ OW am =) < 
eww ee ee een ae. ants a <:-OOa ~N 
pm 8 WE  < COM) Se et tet 4 oO ooa\ovwres 
£5S-250% We 833 4 -wOonrw 
weit octave aoa @rovgm oar N a3 i Ws 
MOW OW SOT VONDONMe OwNoOO w wou OO mm 


VO | 
1 
on 


1. 
| 
| 
| 
| 
| 
| 
| 
| 


) 


ovr 
a 
eo 
a 
os 
a 
= 
ax 
= 
= 
ot 
ox 
‘3 
zo 
oa. 
>u 
_ Soe a=) c 
oo eee —- ODOTVNO NO ec 
te oe | anaacnaoaoaoaa ee OOF OMNI -@eeo 
ev o tvs ) 
oac eere terete tartare rs fe rere terete terete rarare ~w Ow 
=n wv | 
° o os 38 
wu c “avvoo 
= - 8 eo o~ 
>a oe oaeww | 
~~ = «oe “-O | 
<—_—- O&O —- dor Tou —-COOoO MnueAmMmaooewaiwcto N °o -~ oa 
<= SOS o ooa oot MOvrawuSd OOo <a OODvradwwuuwd ~~ — pA | 
> CO u OOou BVuUUKiuuwuod A a OO o +> -- OO 
} oS ms wwe PMMMM—mimriprna MPOMMMMMMMM— MII oO w § Huey 
o o www NOWOAMNMAOD a<cunogaortooaro o =z + @ 
—- O&O o wnw we COOOCMu Oo DLV VVO@ONGAu wo o aocv 
mu Oo oO e&ew woadce-wrectero MOTMVEMoODOODwuO o # u— Owe 
oe «& So eR Y MWOuModor.~TO VCWWOOuwrNMoWwOoo ao « eee & w 
ono oo So wwME “uaQanrOur. fro wiNOO £00 DODu © oS = 5 = 
om S&S oe Whe NOmctODuud OMMMMN aw wh WO So o zuws 
oe OS oe wood VMWOIMNN ee — OO Wr — ANMNoMWwOOauuwo So c¢ -o c 
Ne OS o fer kh —O@Ouwowrnmo WOWNMN a — uu DODWAWNO Go = cc @ - °o 
o- oO oO “Mure Ost wWe WWI OO DVGCLNVOME DMoOwo o -@G8ocvov - #1 
oS oe cvs ww ods umMmMo OA Ousk- DOK. WOMOO oS @ — c  - —J 
vwvTr Oo oS wom VT WINOMDO NMOOwWVOLOOwr tomo oe @ eouewewe® © 3 
ona © oS Ou. VA MDW LvwWO Ve ONnNDWO de ewweo o_— ocwotu - wt ©& 
oo S&S So fOow BNE FMOO~TO OuuroOovrrUMOr. Oo Co = mt 3cc @i ° 
-—- Oo o FOF. DWDLF ACMIDOO MOOT aAOUOO So o ‘ - 30 «e Cc . oO 
ss So oe -Oovw OK LATOR. OS wououddoe IMO o =t weeny QB 
aa & oS wood a AO. OO 7 O NF OMMN OR — hu Wa O oS zetwrm Cc 
ww © o «wom Ooawdwarrr awo ~~ oWOVOOWO Gr. Oo o i ww o 
mm oOo oe oor OOOO wor uO OO-Ow-e- <u ove So c -@o a 
‘6 o oS aot MNDOuUru~TOO Or ru NOr. Od MO o mm @euLl tcc cad 
oo oS So FOU WOVUNM u LO WN fe O- UML OO So « ceo- @ 
- o oe oow ~~ oOouwud-ao <ewe LAO MMvVOO So ‘ ited 
So oe «o- wT Oro oO est OF OMMODWMO o =z oe eo ar wy 
ao oS aoa moor <r-wTwO —uUsn £0004 OM Oro o =z Z2zwtte*-t 8 ££ 1 
o oe «Oct MWMrwNuOr LOO WMNANOVOSCONOVOO o 3 So SBOeMouv ce 
So oS «or MMOuwuwwWusTowo OEEIOMNSWMOUMOOOOWO SO O — # B24 @ 
~N So So mow COWOOaDL OVNOC- WwFTOOovVvOrweOoOmMWwo-—- So a _ otweo = ci 
So So NO MO Qmucw <cO™ We OOF WMOWorr-o~ SO & a mam t«& oOo oF 
o So aOow UNNTH-u ZCOMOA NN &P-DOUVTUICUOOOA OC WwW — «= egos ar o 
=< & =< =< MR >< >< CK OK OK OK OK OK OK ED <> EO OOK KK a wo co ~ x1 
c € cs aac eaceeeece eee e<eeeeeee eee ee Ee « oO vw zcewveo — oi A 
co - lon Yz=ee vw 3S 
cos oOo a w Beewrod wv 
3 i=] -* ail oO O@g@ette © 
o 6 <« «<a q<aq<addeeacaeet i <«<«q4€44¢e ee eeaeaeaw < - aovoa . 
. wee —- eee [all sel eel aan aul eel aed wd aed aed aed aud aed ae eel ed ed ed ed ee DO —_+ *@#—-a @1 
vw Vv VLuU VWVVOVYOUYY + VVVOVVYVUVVOYUYYU + YL WD <«<« »+oax @ vi oOo 
ce oO o ooo oooooooo0ov ooooocoocoeocoo0o0oVv co 2 - ws<t £ @! 
e - eeoneaenen as. Be seeaeen nee aaa . ° Oo ottcu>er - Di " 
— ~N u“ 1) —_ 
Oo — oT oe co 
a - @ z 2 , vw 
= <= w w 4 
oO q — onl _ + p= | 
- Co hn J oO a qa“, a 
ee =z -“—-~N 6 6«O oO =< x = + 
z oa wuw _ —_ w ws ”“” fete terete *e*e tetera teteterterte 
w OM OMPNT NOR. DOO CUM TN OR. DOO K— CIM TN OR. DOHOOK— CIM TWN OR. DOO KK CIPI TN OR. DO OK UIM TNO 
a IIOP SF SF ST SE EP ET TE TNNAAAUAUNOODOOOOO 
RISO OOOO IOI OOO OI OI OI OSI OI OI OI OI OS OS OI OI OS OS OI OI OI OI OI OS OS OI OS OI OI OI OI OI OSI OI OI OI OI ONION OI OI ONIONIONIONS 
ow 
*s 
C= DOOOCOCOCOCOCOCOOOOCOCOOOCOCCOCOOGOOOCOCOOCOCOOOOSOOCOCOOCSCOOCOCOCOOOOOOOOOoOOO 
HK — DLCVVOOwW4 COOK UMTMNAOR.DOOOAS FOV QWKLOK- TUM TNOOORAPPAALALALAL ALLA. 
SSF FF FTF TF TOM AMMAMAMMMAMMMAMIMMAANANUNUANOODOODOODOOOOOOOOOOoOOOOoOOoOoOoOoOooO 
#S OOOO OCOOCOOCOOOOOOOOOOOOSCOOOSCSCOCOOSOSOSCOSCSOSOSSOSOSOCOOSOCOOCCOoOOoOoooo 
wo 
aw oor TMU - OOOe NuAmaocooveiweioa 
w ooa Oo ae OvrOwuOooO Oo <a OOvraoawuu woo 
| i=) oouw. i Ot OO Ae Oe. Oe. Oe Oe. Oe, he, te te te te tc OOO 
vrwre PM PIMMMMmMiImMOoOO PEM MMMMMIMIENOO 
www NOWONIMMASCSO oa<unogoaor coovroo 
monw T- OOOO Mu OO RSSANe Rees stss 
tk et woo-watroO Se Oe Qaucoeo 
Lt ek MOuMaoor.~TrOO Vwwowousy YVwWoooe 
VN uOr-Our. ff <O NCO £00 MODu Oo 
Wn NOMELOOuu OMNAMIMIN uo Ph oo WI 
wood OM WAIN ee Mr — ONONWOOQuu oO 
~~ -O@Muvo wr DOoUNMNua eu OOAWO 
wwe ee ee ee, DVOGVUNVOM cmowo 
wns wwTodtery O Ow. woor. rw a 
oon a Osw Dow NOODwWeDOOyT fUM 
QOouw Ve @Oododcu Vr ONMODO Ee ewe 


row AE FUN OuwuwvOOvrh OMOr. 
or S38 NOOKS Soe Sek 
ov Wu Oud0- = 
Sa S&ananre KFQWOUD OWS 
OOO Wor. OO Owe fu Or. 


Or <r unorOdcerm 


oOo 
AnS<u-S-omaesS 
<hwe LOOK MMvO 
~wTVNST OR on 
WE OD $8 ~ 
VLOoorw 


00000000 00000000 00000000 00000000 
00000000 00000000 00000000 00003F91 
00000000 00000000 00000000 80004072 


CIG2ADIE FA 


SoSSo 
wae 
wh COMUS 


J 2 
= - - e +4 = T 
At ali OTSSPOWMN_RS Hef Loating. #e Heft OSEp-19R4 41:28:21 ATARTL -SREDOTSPORH.MAR;1 °° (hy 1 


POWHH_R3 = H_floating ** H_floating 


28 «656 Undetined Exponent iat ion | 
<0 any Undefined Exponentiation | 
? > ye * Log2(base)) 
> < 2“(exp * log2(base)) 


Floating Overflow can occur. 
Floating Underflow can occur. 
Undefined Exponentiation can occur if: 
ase is 0 and exponent is 0 or negative 
2) base is negative 


The basic approach to computing x**ty as 2“Cy*log2(x)] is the following: 


Step 1: Compute togets) to sufficient precision to guarantee an 
accurate final result (see below.) 

Step 2: Compute y*tlog2(x) to at least the accuracy that log2(x) 
was esnewtae. 

Step 3: Evaluate 2“Cy*log2(x)J] accurate to the precision of the 
datatype in question. 


To determine the accuracy to which log2(x) must be computed to, write 
y*logesx? as I + h, where I is the integer closest to ytlog2(x), and 
= ytlog2(x) - I (Note that thi =< 1/2.) Then 


2*Cytlog2(x)J = 2°¢(1 + h) = (2*1)#(2*h). 


nce the factor 2“*I can be incorporated into the final result by an integer 
dition to the exponent field, we can assume that the multiplication by 

] incurs no error. Thus the total error in the final result is determined 
how accurately 2“h can be computed. If the final result has p fraction 
ts, we would Like h to have at least p good bits. In fact it would be 

ce if h had a few extra guard bits, say 4. Consequently, we would Like 

Ss. 


Let e be the number of bits allocated to the exponent field of the data type 
in question. If I requires more that e bits to represent, then overflow or 
underflow will occur. Therefore if the product y*log2(x) has e + p + 4 good 
bits, the final result will be accurate. This requires that log2(x) be 
computed to at least p + e + 4 bits. 


Since log2(x) must be computed to more bits of precision than is available 


n the base data type, either the next level of precision or multi-precision 
; arithmetic must be used. We begin by writing 


en+1 
log2(x) = log2(b) + Ps c(2n+1) #2", 


n= 
; Where c(1) = 1, and 2° = (2/ln2)C(z-b)/(z+b)]. Hence 
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log2(x) = log2(b) + 2' + > cl2ne1)#2 | 
nzT~ | 
= log2(b) + z* + plz"). 
Note that if p(z') is computed to P bits, and toge (b) + 2' is computed 
to pte+4 bits and overhangs p(z') by e+4 bits, the required accuracy will 
be achieved. Consequently, the essential tricks, are to pick b such that 
the overhang can be achieved and to compute log2(b) + z' top +e + 4 bits. 
CALLING SEQUENCE: 
power.wh.v = OTSSPOWHH_R3 (base.rh.v, exponent.rh.v) 
IMPLICIT INPUTS: 
none 


OUTPUT PARAMETERS: 
none 

IMPLICIT OUTPUTS: 
none 


FUNCTIONAL VALUE: 


The H_ floating result is returned in registers RO-R3. This 
is a violation of the VAX procedure calling standard but is 
allowed for compiled code support routines. 


SIDE EFFECTS: 


Modifies registers RO-R3! 
MTHSK_FLOOVEMAT = Floating overflow , 
MTHSK_FLOUNDMAT = Floating underflow if FU bit is set 
SS$_ROPRAND - Reserved operand fault 
SIGNALS MTH$_UNDEXP (82 = * UNDEFINED EXPONENTIATION') if 
1) base is 0 and exponent is 0 or negative 
2) base is negative 


DDD DDD DDD. DDD. DDD DDD DAP PAA AAA AAA AAA AAA A AAA 
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03F0 $7 0 -ENTRY OTSSPOWHHR3, ACMASK 
67 72 3 
67 6 ; Move x to RO/R3. If x < 0, or x = 0 and y =< 0, return ‘UNDEFINED 
of iS : EXPONENTIATION’ error condition, otherwise attempt to compute x*tty 
67 ik: : 
E 34 ¢2 067 7 SUBL #52, SP ; Allocate space on the stack 
50 "04 AC 70FD 0675 378 MOVH = base (AP), RO t RO/RS <a= x — 
1D 14 Q67A 7 BGTR DEF INED : If x > 0 attempt to compute x*tty 
07 19 067C 0 BLSS UNDEF INED ; Branch to error code for x < 0 
14 AC 73FD 43 1 TSTH te Bad. 3; Check sign of y (Mote that x = 0) 
01 15 debe BS BLEQ UNDEF INED ; Branch to error condition if y =< 9 
684 84 ; 
Dee $3 ; If processing continues here, this implies that x = 0 and y > 0. Return 
0684 ; xeey = 0 
gi: y 
oes, | 368 
04 0684 89 RET ; Return 
bess 391 
0685 3 : If brorese tne continues here, this implies that an undefined exponentiation 
0685 93 ; was attempted. Signal error and return 
Bs 
bege 38 UNDEF INED: 
50 7CFD 0685 9 CLRO RO H 
50 000 8F BO 0688 98 MOVW #*x8000, RO 3; RO/R3 <-- Reserved operand 
7E OO°'BF 9A 0680 99 MOVZBL #MTHSK_UNDEXP, -(SP) ; Put error code on stack 
00000000'GF 01 FB 0691 400 CALLS #1, G*MTHSSSIGNAL ; Convert error number to to 32 bit 
0698 401 3 condition code and signal error. 
0698 106 : NOTE: Second argument is not re- 
boas 40 : _ quired since there is no JSB entry. 
04 069 404 RET ; Return 
$699 406 
0699 407: If proqssscng continues here will attempt to compute x*ty as cory toge tnd. 
0699 408 ; We begin by determining k and f such that x = 2“k*f, where 1 =< f < 2. 
0699410 
0699 411 DEFINED: 
54 50 FFFF8000 8F CB 0699 tig BICL3 #*XFFFFB8O00, RO, R4 : R4 <-- biased exponent of x 
54 00004001 Hy cs 6A1 41 SUBL #*X4001, R4 3 R4 <-- k = exponent of x - 1 
50 4 ¢C 6AB 414 SUBL R4, R 3 RO <-- f = 2e(fraction field of x) 
gas ae 
6AB 119 : We are now reedy to compute log2(x). This computation is based on the 
ny $18 3 following identity: 
6AB «420 ; eerie. ' f-a 
6AB 421; Log2(2*k*f) = k + log2(a) + ----- > won- 2°(2j+1), where 2 = ----- ‘ 
m4 ? ¢ : in(2)/____ 2jet fra 
6AB 424 : We begin by determining a as b*i, where b = 2°(1/64) and i is 0, 2, 4, ... 
6AB 425 : 64 or 1. pecifically i is chosen by table look-up in such a fashion as 
6AB 426 ; to minimize the magnitude of z. Since log2(a) = 1/64 we may write 
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AB 427 ; 
. ? 8 3 log2(x) = k + 1/64 + z2*p(z*z). 
AB 450 ; NOTE: For i = 2, 4, ..., 64, we may write i = 2n, and hence 1/64 = n/32, i.e. 
AB 431 ; a is an integral power or 2°(1/32). These values are stored in Al TABLE and 
AB 4 § : AG TABLE. For i = the value of 2*(1/64) is stored janed iegely BEF ORE 
6AB 435 ; AI_TABLE and A2_TABLE. Consequently, to access the value of 2*(1/64) from 
6AB 434 ; the table, a negative index is used. 
gaa a3G 
6AB 4 ; EVAL_LOG2: 
3s 30 7 % 6AB 4 3 ROTL #7, RO, RS : R5(0:6) <== high 7 fraction bits of f 
54 554 ea. 6AF 4 ASHL #6, R4, R4 3 R4 <== 2°6ek 
55 FFFFFF8O 8F CA 06B3 440 BICL AXFFFFFFB8O, RS ; RS <== index to INDEX table 
55 FFFFFOSF ore? 3 ra 723 syst g MOEKERS I, R5 3 oe = 1 of -1 
: Branc 
9° D6 06C4 108 INCL R4 7 R4 <m= 2°6(k + 1/64) 
55 D7 06C6 444 DECL R5 3 RS <-= =2 
03 11 Q06C8 445 BRB 2$ ; Join common code 
54 55 CO Q6CA 446 1$ ADDL R5, R4 3 RS <m= 2°68(k + 1/64) 
Oecd a8 
06CD 449 ; Since there is no back up data type to compute the necessary goers bits, we 
06CD 450 ; proceed by computing 2 = (f-a)/(f+a) in t¥e pores: z= 21+ 22, where z1 is 
06CD 451 ; the high 53 bits of z and 22 is the low 11 its of z. Further, to obtain 
06CD 12g ; the desired accuracy it is_necessary to work with a = al + a2, where al and 
06CD 453 ; a2 are the high and low 113 bits respectively of a. We begin 
Beco tee 3 computing (in G-format) 
Beco $28 3 z1 = (f = al)/(f + al) 
06CD 458 ; Note that f-al can be computed exactly in 113 bits, but f+al may require 114 
06CD 459 ; bits. The 114th bit can be determined by the exclusive or of the low bits of 
06CD 460; f and ail. 
gh i8) | 
55 55 FF 8F 78 0O6CD 163 2s ASHL #-1, R5, RS ; RS <=-= octaword offset into Al, A2_TABLE 
56 FEEFF9RG EF4S 7DFD O6b2 464 MOVO. L*A{_TABLECRSJ, R6 > R6/R9 <== al 
sc 5 59 CD 060DB 465 XORL3 RY, R3, (SP) ; SP -=-> XOR of low bits of al and x 
O6DF 466 : (This will be used to determine 
sed 467 3 the 114th bit of f+al.) 
04 AE 50 56 63FD O6DF 468 SUBH R6, RO, t2(SP) 3 t2 <-- f = al (exact) 
14 AE 50 56 oiF8 6E5 rt ADDH R6, RO, t4(SP) 3 t4 <-- f + al (rounded) 
26 14 AE 76FD Q6E & ? CVTHG t4&(SP), R 3 R6/R7 <-- f + a | 
5 04 AE 76FD O6F 47 sy ues t2(SP), R 3: R8B/RO <-- f - al 
58 28 46FD O6F t76 DIVG R6, RB 3: RB/RY <== 21 (G6) | 
24 AE 8 56FD O6F 47 CVTGH RB, t6(SP) : t6 <== 21 (H) 
ere 98 
ore 478 ; To compute 22 we note 
re $78 ; z= 21+ 22 = (f = al = a2)/(f + al + a2) 
ore ; 9 : ==> 22 = (f - al = a2)/(f + al + a2) = 21 
6FE 4 : : Now let v = f + al + a2 = vl + v2, where v1 and v2 are the high 49 and low 
6FE 483; 115 bits of v respectively. Then | 
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é : z2 = C(f = al = zi*vl) = (a2 + 21*v2)I]/v 
2 : We begin by computing v1 and f = al = z1*v1 
6 e 
56 4A 7D MOVQ t4(SP), R6 3 R6/R7 <== high quadword of f + a1 
58 06 DOP mova. #0, S/N ae wi 
50 14 AE 6 63FD SUBH R6, t4(SP), RO 3 RO/R3 <-- w- vl (exact) 
56 24 AE 64FD MULH2 t6(SP), + R6/RO <== z1"v1 (exact) 
62FD SUBH R6, t2(SP) 3 t2 <-- f - al = zitvl (exact) 


: Compute v2 and a2 + altv2 
BICL @*XFFRFEFFFF, (SP) 


6E FFFEFFFF 8F Check if w was rounded 


CA : 
13 BEQ ; Branch if not rounded 
50 FD9D CF oer? SUBH TWO_M112, RO 3; Correct for rounding error (exact) 
50 FFFFFB84 EF4S 60FD 3$ ADDH L“A2_TABLECRSJ, RO 3: RO/R3 <-- v 
50 24 Af 64FD MULH t6(SP), 3 RO/RE <== 21*y2 
50 FFFFFB76 EF4S 60FD ADDH L“A2_TABLECR5], RO 3 RO/R3 <== a2 + z1*v2 


: Compute 22 


SUBH RO, t2(S 3 t2 <-- (feal-z1*v1)-(a2-z1*v2) 
DIVH tatSP), st2 (SP), R6 3 RE/R9 <== 22 


ns 
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04 AE 50 62FD 
56 04 AE 14 AE 67FD 


The next step is to compute log2(x) accurate to at least 128 bits. This is 
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7 
7 
7 
7 
7 
7 
7 
07 
07 
Bf 
7 
07 
07 
8? 
7 
07 
0740 
0747 3 
0747 : 
8) 3; accomplished as follows, let 
0747 3 w= Fake tppece x) 
Bf 7 : = (2°6)Ck + 1/64 + gsptzez)3 
747 : = 2°6(k + 1/64) + (2%6)*z"(c0 + c2*2*3 + 1... + €10*2°11) 
0747 : = isa) + 1/64) + 2°) + 2'(d2"z'*2 + ... + d10*z°*10) 
0747 : = C2*6#(k + 1/64) + 2°) + 2z*eQ(z'#2") 
Bf : = wi + we 
0747 3 where z' = (2*6*c0)*z and wl and w2 are. the high 49 and Low "ss bits of w 
0747 3 respectively, Note that the choice of ‘a’ used in computi ng 2. uarantees 
0747 3; that 2° oyerhgngs z'tq(z'*z') by at least 15 bits. Hence, if w is computed 
Boe) 3 as wi + w2, 128 bits of accuracy can be obtained. 
15H 3; We begin by defining 
747 ; c. = high 113 bits 9 of 5h 0e8?? 
747 : c1 = high 49 b (2*6*c0) 
747 3 c2 = low 113 bits _ (2*6*c0) 
747 3 then 
747 ; z* = (21 + 22)#(cl + ¢2) 
747 ; = z1*cl + 21*c2 + z2tc, 
ee 
50 26 AE FDA4 CF 65FD 0747 MULH3 C2, t6(SP), RO 3 RO/R3 <== c2*z21 
24 AE FDC CF 64FD 74F MULH C1, t6(SP) 3 t6 <== cl*z1 
56 F075 fF 64FD 756 MULH C, R6 3 R6/R7 <== #22 
AE 56 50 61FD 075 ADDHS RO, Rb, t4(SP) i th <om cazd + c2ezi 
56 14 AE 24 AE 61FD i ADDHS t6(SP), t4(SP), R6 > RO/R9 <== 2° 


soem 


52 26 


wu 

No 
MACUL SNS 
COMO WwW Om 
a Ie I ent en Be Bene Bo 
cvoVvVvVT7TTC9o 


2 
14 AE 24 
04 AE 
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14 AE 5 


70F D 


54 50 


FFFF8000 8F 
4000 8F 
55 04 AE 


C 
A 
FFFF8000 BF C 
55 4000 8F : 
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55 FFF 8F Bi 
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We proceed by letting 
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wi = high 49 bits of 2*%6*(k + 1/32) + 21*¢1 
w2' = {C2*6e(k + 1/32) + z1*c¢1 = wi) + 21*¢2) + 22ec,. 


and 


==> 2*6e(k + 1/64) + 2' = 


==> w 


wl + we, 


where w2 = w2" + 2°#Q(z'*z') 


CVTLG RG, RG 
CVTHG t6(SP), R2 
ADDG3 R4, R2. RO 
SUBG3 R4. RO, R2 
CVTGH ‘ 
SUBH2 R2. t6(SP) 
ADDH2 té6(SP), t4(SP) 
CVTGH RO, t2(SP) 
Compute w2 
MULH3 = R6, -R6, RO 
POLYH RO. #LOGLEN-1, LOGTAB 
MULH2 R6, R 
ADDH2 RO. t4(SP) 


We now calculate yetogets) = Cyleyg 
yl and y2 are the high 56 and 


MOVH exp(AP), RO 
reat tor the possibility of overflow in the computation of 
S Ww 
than 16383. 


3 #*XFFFF8O00. RO, R4 


BICL 
SUBW2 R4 
BICL a-xt FFr6000, t2(SP), 
SUBW2  #*X4000, R 
ADDW2 = R4, 
CMPW OW“ K3FFF, RS 
BGEQ NO SYS vER FLOW 
BRw Y_TIMES_wWi_OVER 

95 NO_SYS_OVER FLOW: 

36 MULA2 RO, t4(SP) 

9 MOV RO. R4 


wl + 


e(wl )= 
ow 57 bits of y 


l_occur if the exponent of y plus the exponent of wi 


R5 


we’. 


C2*6e(k + 1/64) + 2°) + 2z*8Q(z'*2") 
wi # w2' + 2'#q(z'#2') 


R4/R5 <-- 
R2/R3 <-- 
RO/R1 <== 
R2/R3 <-- 


2*6(k + 1/64) 
z1*c 
2*6(k+i/32) + 21%¢1 


R2/R5 
C2*6(k+i/32) -wl+z1*¢2) 
t4 <== w2' 
t2 <== wl 


RO/R3 <== 2'*2' 

RO/R3 <== q(z'*z') 
RO/R3 <-= 2°2Q(2z'*2') 
t4 <== w2 


+w2 yitwl + y2ewl + yew2, where 


respectively. 


; RO/R3 <-- y 


twi. 
is greater 


biased exp of 
unbiased exp of y 
biased exp of 
unbiased exp o 
unbiased exp of wit 


largest unbiased ome possible is 16383 


4 <-- 


; % y*we- . 
>; R4&/RS <== high 49 bits of y 


Page ig, 


bits of z1*c1 included in wi 
<-- Ne of 21*c1 included in wi 


. 2 
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56 52 FFFFOIFF 8F cB 7DA 238 BICL3 #*XFFFFOIFF, R2, R6 3 R4/R6 <== high 56 bits of y 
7 00 OD rEg 9 MOVL. #0, R? : RG/R? <e= yl | 
0 4 62FD 7E 600 SUBH R4, RO 3 RO/R3 <-- yg 
50 04 AE 64FD O7E9 601 MULH2 t2(SP), RO > RO/R3 <-= yoew! 
04 AE 54 64FD O7EE 6 é MULH2 RG, t2{SP) t t2 <e= ylewl 
oe 0a | 
er oR? 3; The next step in computing 2“Cy*log2(x)] is to write y*log2(x) as 
Aa $07 3 ytlog2(x) = 1 + j/32 + g/32, 
7F 609 ; where I is an integer, j is an inte s" petyeen 0 and 31 inclusive, and 
7F 610 ; g is a fraction in the interval [-172, 1/2) 
OrFs eID | 
54 14 AE 04 AE 61FD grr gi8 ADDH3 t2(SP), t4(SP), R4 3 R4&/R7 <m= ylewl + yew? 
4 43h CF 60FD 7FA 6146 ADOH SHIFT, R4 : 
54 FE5B gf rere 0800 615 SUBH SHIFT, R4 3 RG&/R7 <== 2°6(1 + j/32) 
04 AE 4 62FD 0806 616 SUBH RG, te(SP) 3 t2 <-- those bits of z1*y1 not 
080B 617 : included in 2°6(I + j/32) 
50 04 AE 60FD O80B 618 ADDH t2<SP), RO : 
50 14 AE 60FD pety 619 ADDH t4(SP), RO : RO/R3 <== 2*7(g/32) } 
58 54 8000 8F AB 081 620 BICW #*x8000, R4, RB ; RB <== exponent field of 2°6(1+j/32) 
54 D7? O81B 621 DECL RG : R4 <m- 2*52ytlog2(x) 
58 4013 8F B1 081D 6 ¢ CMPW #*x4013, RB ; 
36 15 0822 «66 BLEQ EXCEPTION_1 : 
58 54 6AFD 83 3 ° ¢ CVTHL R4, RB 3 RB <-= 2°5%(] + j/32) in integer 
0828 626 ; 
oe ; 6 4 3; We can now compute 
0838 ? ; : xtty = 2*Cy*log2(x)] = 2°CI + j/32 + g/32) | 
0898 o3! $ = (2°1)*CA(B+1)] = 2*I1*CA + A*B), where 
0828 $38 : A = 2°(j/32) is obtained by table look-up and B = 2*(g/32) - 1 is obtained 
0828 $3 3; by a Min/Max approximation. 
0858 636 | 
FD61 CF oc 50 75FD 0828 637 POLYH RO, #EXPLEN-1, EXPTAB 3 RO/R3 <-- B = Ss bef - 1 
59 = 558 FFFFFFEO 8F CB O082F 638 BICL3 #*XFFFFFFEO, R 3; RY <== index into Al_TABLE | 
59 F853 CF49 7EFD 08 7 639 MOVAO Al TABLELRO, R9 : RO <== address of A | 
50 69 64FD 083 640 MULH (RJ), RO : RO/R3 <== AtB 
50 9220 C9 60FD 084 641 ADOH TABLEN(R9) , RO : RO/R3 <== Ase + A2 | 
69 60FD 4 ot¢ ADDH (RI), RO $ he <-- 2*C(j+g)/32] = (A*B+A2)+A1 
58 1F) =O AA (BKC OG BICw #*xiF, RB : R7 = 2*5e] 
58 58 FB ef 9C O84F 644 ROTL #-5, RB, RB : | 
50 8 *9 54 645 ADDW R8, RO 3 RO/R3 <-= 2*1*#2*C(j+g)/32) 
o8 1 57 966 BLEQ EXCEPTION 2 ; see what exception is if neg or = 0 
04 59s «64 RET ; otherwise return result in RO 
ett: 
* 620 ; Handlers for software detected over/underflow conditions follow | 
5A 6 § 
SA 653 EXCEPTION 1: ightep 
50 73FD 5A 654 TSTH RO : if big ARG > 0 goto overflow 


dD 3 
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6 ee BGEQ OVER 3 handler, otherwise go to 
09. BRB UNDER ; underflow handler 

EXCEPTION 2: 
57 65 TSTw R7 ; test sign of I; if I <0 
19 18 BGEQ OVER ; go to overflow handler 


; y*wl would have caused a hardware system floating oyer tion error. If y<0, 
; then we should return a result of 0 since result = 2*(y*(wltw2)). Note, 
; y can not be zero. 


Y 


_TIMES_W1_OVER: 
50 73FD TSTH RO : if y < 0 no overflow is needed 
14 16 BGTR OVER : overflow | 
UNDER: 
50 7CFD CLRO =. RO ; RO/R3 <== 0 
08 04 AD 06 «€E1 BBC #6, SFSW_SAVE_PSW(FP), 2$ 


has user enabled floating underflow? 


7E O0°SF 9A MOVZBL M#MTHSK_FLOUNDMAT, -(SP) : Put underflow code on s 


tack 
convert to MTH$_FLOUNDMAT (32-bit VAX-11 


Q0000000'GF 01 FB 


; exception code) 
: signal condition 


i #1, G*MTHS$SIGNAL 


C 
2$: R 


0 

4 

9 

’ 

§ : Underflow; if user has FU set, signal error. Always return 0.0 
4 

8 

9 

0 

1 

‘ return 

3 

8 : Signal floating overflow, return reserved operand, -0.0 
8 


7E OO°'SF _ 9A 
50 8000 8F B60 


oO 
~N 
Oo 
pms 
o 


R RO = result = reserved operand 
MOVW #*x8000, RO “0.0. RO will be copied to 
signal mechanism vector (CHFSL_MCH_RO/R1) 
so can be fixed up by any error 
handler 
signal condition, 
return - RO restored from CHFSL_MCH_RO/R1 


00000000°GF 01 FB 


695 cms #1, G*MTHS$SIGNAL 


OOO O000 90 9 NI NI SNS SI SINS SSN NPP PPD AP AA AAA AAA AA AAA 
MRM — PP YS FS VIIVUMMMMMMOeeXINMS PY} YS \S PY YF OMIM 


SOOCCOCOCOOOOOOOOOOCOCOOCOCOOCOOOOOOOCOOOOCOoOO 


9 OVER: ot ee #MTHSK_FLOOVEMAT, -(SP) ; Put overflow code on stack 
698 ~END | 
| 


OTSSPOWHH 
Symbol table 


A1_TABLE 
A2- TABLE 
ACMASK 
BASE 


OvE 

SF $W_SAVE_PSW 
SHIF 

14 

16 

TABLEN 
TWO_M112 
UNDEF INED 


UNDER 
Y_TIMES_W1_OVER 


PSECT name 


Initialization 
Command processing 


Pass 

Symbol table sort 
Pass 

Symbol table output 


ite Biabne 


Page faults 


~ REAL*16 ** REAL*16 power routes : i ak i Sat 9 96:00:37 


90 R 
mig 
= F 
= 04 
4D0 R 2 
4EO R 
4FO R 
699 R 
sete 
SAR 0 
seis e 
= 
= $88Rh 
00000 3 R 88 
00000000 R 0 
= a4 4 4 
00000500 R 02 
eeeeeene x 0 
eeeekeee x + 
eeereeee x 00 
keekeeee x 00 
000007D2 R 02 
§0090676 RG 02 
0000087E R 02 
= 00000004 
00000660 02 
= 00000004 
= 00000014 
= 44 ts 
= 00000220 
000004C0 R 02 
00000685 R 02 
0000086A R 02 
00000865 R 02 
teeeoeeeoeeeeoeeces ea 
! Psect synopsis H 
Allocation PSECT No. Attributes 
00000000 <( 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE 
00000000 ( 0.) O10¢ +1.) NOPIC USR CON ABS LCL NOSHR- EXE 
00000892 ( 2194.) 02 ¢ 2.) PIC USR CON REL LCL SHR EXE 


$woceccceceeeeeoeocececen} 


! Performance indicators : 
pew oeeeeeoeeeeeeeceseces 


CPU Time 


8:00: 2:96 
HY 
8 


0. 


AX/VMS Macro 
MTHRTL.SRCJOTSPOWHH.MAR; 1 


v04-00 


NORD NOWRT NOVEC BYTE 
RD WRT 
RD NOWRT 


NOVEC BYTE 
NOVEC QUAD 


wef 
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Psect synopsis output :00: +e :00: Be 
Cross-reference output :00:00.00 :00:00.00 
Assembler run totals 43 :00:05.16 :00:17.64 
The working 8 Limit was 1200 pages. 

13981 bytes (28 pages) of virtual memory were used to buffer the intermediate code. 
There were 10 pages of symbol table space on ecares to hold 62 non-local and 4 local symbols. 


758 source Lines were read in Pass 1, producin object records in Pass 2. 
pages of virtual memory were used to define & macros. 


$eeroceeeec oe ececea Sees es ecaas + 

! 

$eemmme wwe wwe eee eeweececee= + 
Macro Library name Macros defined 


-$255$DUA28: CSYSLIBISTARLET.MLB; 2 4 
88 GETS were required to define 4 macros. 


There were no errors, warnings or information messages. | 
MACRO/ENABLE=SUPPRESSION/D1I SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWHH/OBJ=0BJ$:O0TSPOWHH MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MSRC 


| 0265 AH-BT13A-SE ) NT CORPORATION 
gave VAX/VMS V4.0 D PROPRIETARY 


